Multi-orbital physics in Fermi liquids prone to magnetic order 
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The interplay of spin-orbit-coupling and strong electronic correlations is studied for the single-layer 
and the bilayer compound of the strontium ruthenate Ruddlesden-Popper series by a combination 
of first-principles band-structure theory with mean-field rotationally invariant slave bosons. At 
equilibrium strongly renormalized (spin-orbit-split) quasiparticle bands are traced and a thorough 
description of the low-energy regime for the nearly ferromagnetic bilayer system in accordance with 
experimental data is presented. The metamagnetic response of Sr3Ru207 in finite magnetic field H is 
verified and a detailed analysis of the underlying correlated electronic structure provided. Intriguing 
multi-orbital physics on both local and itinerant level, such as e.g. competing paramagnetic and 
diamagnetic contributions, is observed with important differences depending on the magnetic-field 
angle 6 with the crystallographic c axis. 

PACS numbers: 71.27.+a, 71.18.+y, 75.10.Lp, 75.70.Tj, 75.47.Lx 



I. INTRODUCTION 

Ternary ruthenium-based oxide compounds condense 
in a rich variety of different crystal-structure types with 
associated rather delicate electronic and magnetic prop- 
ertiesii Within the perovskite-like Ruddlesden-Popper 
series of the ruthenates A„4.iRu„03„+i (A=Sr,Ca), 
where n labels the layers of corner-sharing RuOg octa- 
hedra (separated by SrO or CaO rocksalt layers), the 
interacting electrons pose a specifically challenging prob- 
lem. In that family of seemingly rather similar com- 
pounds, intriguing competitions between Fermi-liquid, 
Mott-insulating and superconducting behavior occur in 
conjunction with particular complex magnetic response. 
For instance, albeit Ca and Sr are isovalent, the re- 
spective single-layer compounds A2RUO4 exhibit dras- 
tic different phenomenology, since Sr2Ru04 displays un- 
conventional superconductivity below Tc^l.b K^^— , while 
Ca2Ru04 becomes an antiferromagnetic (AFM) Mott in- 
sulator below Tmi~77 Ki^i^ As an important generic dis- 
tinction, compared to the Sr compounds the Ca subclass 
in the series exhibits stronger distortions and deviations 
from an ideal high crystal symmetry. 

A prominent aspect of the intricate physics in the over- 
all metallic strontium ruthenates is the onset of ferromag- 
netism with n. Concerning the series end members, the 
perovskite SrRuOs (n— >cxj) is ferromagnetic (FM) below 
Tc~165 K and tetragonal Sr2Ru04 ('^.= 1) is paramag- 
netic (PM) at ambient temperature T, but shows FM 
tendencies as well as incommensurate spin fluctuations at 
q=(±0.67r/a,±0.67r/a,0)^^ The (n=3) Sr4Ru30io and 
the (?T-=2) Sr3Ru207 compounds are both orthorhom- 
bic, but whereas the former is verified FM)^ the latter 
is still PM down to low temperatures. However the bi- 
layer system appears to be located rather close to the 
transition towards FM order^iii with puzzling metamag- 
netic (MM) behavior in applied field below Tmm~ 1 K 
(see Ref. |Tl| for a recent review). As revealed from de 
Haas- van Alphen (dHvA),— angle-resolved photoemis- 
sion spectroscopy (ARPES),— "— opticsi^ and resistiv- 
ityiS measurements, the n=l,2 compounds (see Fig. [1]) 



both belong to cases of quasi-twodimensional (2D) elec- 
tron systems, i.e., show a strong anisotropy between 
transport in the a6-plane and along the c-axis of the crys- 
tal structure. Signatures of strong electronic correlations 
are compelling for these layered ruthenates, e.g. from 
large mass renormalizationS )^°'^^i^^ and originate from 
the less-screened Coulomb interactions within the Ru(4d) 
shell. Nominally four electrons occupy this 1=2 manifold, 
i.e. an Ru'*+ oxidation state may be assumed. The meta- 
inagnetism with applied field H in Sr3Ru2 07 is well docu- 
mented by a large slope ^^f I/^mm around Hmm=5-5{7 .7)T 
for H\ |a6(c)i^i2i Furthermore this MM region may be as- 
sociated with being in the neighborhood of a quantum- 
critial point that can be approached via tuning the po- 
lar angle 9 between magnetic field and the c-axis.— 1^ 
Reachable within fields H<10 T, the MM phenomena 
are acting on a very low energy scale of the order of at 
most a few meV. In this respect, the Fermi-liquid regime 
in vanishing field exists below 10-15 K, however can be 
driven to zero temperature with applied fieldi^ 

In this work we present a theoretical investigation 
of the n=l,2 compounds in the normal state of the 
low-temperature regime with an emphasis on the in- 
triguing physics of the bilayer system in applied mag- 
netic field. The study is based on the combination 
of a first-principles band-structure approach, spin-orbit- 
interaction treatment in the Russell-Saunders limit and 
mean-field many-body theory. The peculiar low-energy 
physics of the layered ruthenates ask for a very detailed 
examination of the single- and many-particle terms in 
the Hamiltonian in order to capture the important pro- 
cesses that drive the physics of these systems i^i^ d^i^^" — 
Besides the thorough description of the crystal bond- 
ing which leads to a multi-orbital based band manifold 
at the Fermi level, it was show n^^i'^°i^^ that addition- 
ally spin-orbit effects play a vital role in the low-energy 
regime. Focussing on the many-body part, e.g. the rele- 
vance of the Hund's coupling Jh in addition to the larger 
Hubbard U was elucidated in several works j^^J^ii^i Much 
theoretical effort has also been devoted to the descrip- 
tion of the MM phenomena in Sr3Ru207, either based 





FIG. 1. (Color online) Crystal structures of Sr2Ru04 (left) 
and Sr3Ru207 (right). Large grey: Sr, blue (dark): Ru and 
small red (dark): O. 



on effective single-band modelings^^""— i^ or with includ- 
ing multi-band degrees of freedom.— "—i^ The existing 
model studies are able to account for the principle ap- 
pearance of metamagnetism, often accompanied by ne- 
matic orderj^l i.e. broken rotational symmetry. Impor- 
tant ingredients for the MM behavior are van-Hove sin- 
gularities close to the Fermi surface, already revealed in 
the single-band approaches!^ More sophisticated multi- 
orbital investigations have been employed to discriminate 
between the importance of dxzi dyz (formally quasi-lD 
like dispersions) and dxy (formally quasi-2D like disper- 
sions) orbital degrees of freedom together with spin-orbit 
coupling (SOC). However often the broken fourfold sym- 
metry in Sr3Ru207 is neglected, when modeling the elec- 
tronic states. In addition, most theoretical works treat 
the many-body interactions in the Hartree-Fock approxi- 
mation, not allowing for explicit self-consistent renormal- 
izations and ill-defined for the metallic regime. 

A combination of the local density approximation 
(LDA) to density functional theory (DFT) with slave- 
boson^"— theory in a rotationally invariant representa- 
tional^ was utilized in Ref. |4a in order to account for 
renormalized quasiparticle (QP) behavior. Here an ex- 
tension of that work is provided by including the effect of 
the spin-orbit coupling on the Fermi-liquid regime of the 
layered ruthenates including applied magnetic fields. In 
the case of Sr3Ru207 the main focus is on a many-body 
modeling that starts from the realistic low-symmetry 
band structure. Because of the fact that the MM prob- 
lem involves very low-energy scales in an underlying low- 
symnretry lattice, we believe that the bilayer system 
serves as a challenging test case for the reliability and 
accuracy of current extended-LDA approaches. 



II. THEORETICAL APPROACH 

The LDA part of this work is performed using an im- 
plementation^ of the highly-accurate mixed-basis pseu- 
dopotential (MBPP) technique^^ employing normcon- 
serving pseudopotentials^ and an efficient combined ba- 
sis consisting of plane waves and additional localized or- 
bitals. For the inclusion of spin-orbit coupliirg and to 
treat local many-body interactions, the realistic Hamil- 
tonian 
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is used, where k denotes the wave vector, a numbers the 
unit cells with i,j marking the Ru ions within, the spin- 
projection is given by cr=t, J, and rf'^' annihilates (cre- 
ates) electrons in the t2g-like Wannier orbitals m, m' . 
The Kohn-Sham dispersion £k in the latter basis is here 
obtained from a maximally-localized Wannier-function 
(MLWF) construction^2iS based on the band structure 
revealed from the MBPP calculation. The prime in ([T]) 
indicates that strictly local (on-site) contributions are ex- 
cluded in the fc-dependent dispersion (see eq. ([S])). More 
details on the determination of the dispersive part in the 
case of Sr3Ru207, where the band Hamiltonian amounts 
to an 12x12 matrix due to the fact that the primitve 
unit cell encloses four Ru ions, is provided in Ref. |4^. 
For Sr2Ru04 the primitive unit cell contains only one 
Ru ion and the Kohn-Sham problem asks for the diag- 
onalization of a 3 x 3 matrix. The unit-cell Hamiltonian 
T-La'^ decomposes into the four terms, reading 

^i'°'' = ^i"'^ + ^i"'^ + ^i""' + ^2"'^ (2) 

and is evalutated in each particle sector of an effective 
i2g problem. The first contribution includes the on-site 
crystal-field through 
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with e\J^, =1/A^k Ek ^Limm' computed from the com- 
plete Kohn-Sham dispersion. We continue with the spin- 
orbit interaction Ha of Russell-Saunders type (or LS 
coupling scheme) on each of the rather light Ru ions. 
Summing over the individual ion contributions leads to 
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where A is the spin-orbit coupling parameter, L^ the to- 
tal orbital momentum operator. S^ the total spin opera- 
tor and ii=\ji-\Si the total angular momentum operator 
for each Ru ion, respectively. Note that in the Russell- 
Saunders approximation these three operators are true 
many-particle operators, given by the sum over the re- 
spective operators for each individual electron p within 



the t2g shell. 



e.g. 
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for the total orbital mo- 



mentum operator. For the t2g orbitals a valid choice 



for the matrix elements {m \Li 
resolved^ 



reads component- model, which reads 
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Due to the cubic crystal-field terms in the ruthenates, 
this representation is obtained by truncating the appro- 
priate matrix elements of a full 4(i shell based on cu- 
bic harmonics to a pure t2g shell including the states 
'idxz^ '^dyz and 4,dxy Although there is some minor inter- 
mixing with the Eg states at low-energy for Sr3Ru207 (see 
section nil)) , that approximation proves to be adequate on 
the present level of the investigation. 

The third contribution to eq. ^ describes the local 
interaction of the Zeeman type with a magnetic field H 
and can be written as 
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Notice that in the weak-coupling regime |H|^A consid- 
ered in our calculations, the spin-orbit coupling domi- 
nates the magnetic-field interaction. Thus eigenvalues of 
Lz and Sz are no good quantum numbers of the system, 
but { J^, Jz} are now commuting with lia . Therefore 
we have to perform the projection of L^ , S^ onto J^ ac- 
cording to the Wigner-Eckart theorem using the opera- 
tors' common eigenspace representation. The term 'Ha '" 
has then the following form: 
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where the notation {■■■)i^sj indicates the expectation 
value in that eigenspace defined by the eigenvalues of 
the operators L^, S"^ and J^, labelled by the quantum 
numbers L, S, and J, respectivelyi^ Notably, the ob- 
ject gi(LSJ) is the generic matrix representation of the 
Lande factor (or g-factor) in that eigenspace. After that 
eigenspace computation one has to transform 'Hq " back 
into the original Hilbert space, where all other parts of 
the local Hamiltonian were derived. Note that in this way 
the g-factor is calculated seperately for each considered 
state. 

Last but not least, eq. ^ includes the electron-electron 
interaction Ha provided by a multi-orbital Hubbard 
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with n=d^d. In eq. ([5]) the first term marks the intra- 
orbital Coulomb interaction with Hubbard U and the 
second term provides inter-orbital Hund's rule corrected 
interaction with t/'=C/— Jh and [/"=[/— 2 Jh for unequal 
and equal spin projections a, respectively. The last two 
parts account for spin-flip and pair-hopping processes, 
vital to enforce the rotational invariance. These terms are 
especially important concerning the magnetic response of 
an interacting systemi^ 

The complete Hamiltonian ([Ij embodies three inter- 
action parameters, namely the spin-orbit interaction A, 
the Hubbard U and the Hund's exchange Jh- For both 
systems the value A=0.09 eV, as obtained from LDA cal- 
culations for Sr2Ru04 by Haverkort et al^, is used for 
the SOC. Concerning the Coulomb interactions, previous 
^Qj,]jg25j26^j3i^ located the Hubbard U for the layered 
ruthenates in the region 1.5—3.1 eV. Here we choose the 
moderate value U=2 eV. With including SOC, that order 
of magnitude is sufhcient to account for the key renor- 
malization effects at low energy. The Hund's exchange is 
fixed to Jh=0.35 eV^i^i^i^ throughout this work. 

Our interacting ruthenate problem is solved via the 
rotationally invariant slave-boson (RISE) formalism^i^ 
in the saddle-point approximation. It amounts to 
a decomposition of the electron's QP (fermionic /j.^) 
and high-energy excitations (taken care of by the set 
of slave bosons {(f>}) on the operator level through 
d^^=-R[{0}]J^J^, /i/'(j', where z/ is a generic orbital/site in- 
dex. Additional constraints for the normalization and to 
match the fermionic and bosonic contents are enforced 
on the mean-field level. The RISE electronic self-energy 
S(a;) at saddle-point is local and incorporates terms lin- 
ear in frequency as well as static renormalizations. It is 
thus given by 



S(a;) =Lu{l-Z- 



with 



h S"*^' , (9) 



whereby Z is the QP-weight matrix and A describes the 
matrix of Lagrange nmltipliers for the enforcement of the 
constraints. Expectation values of any given local oper- 
ator O may be computed via the slave bosons according 
to 



{0)=Y,{A\0\B)Y,' 
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with A, B denoting atomic states and q as the QP in- 
dex. For more details see Ref. |5j. In the present scope 
the method may also be interpreted as a simplified ap- 
proach to solve the dynamical mean-field theory (DMFT) 



equations (see e.g. Ref. lo for a review), compared to 
e.g. more elaborate quantum Monte-Carlo (QMC) tech- 
niques. Since we are interested in the low-energy physics 
of the layered ruthenates at rather small temperatures 
(where QMC usually becomes very challenging), this 
approach is thus well suited to access the Fermi-liquid 
regime including its extension to magnetically ordered 
phases. 



III. CORRELATED ELECTRONIC 
STRUCTURE OF THE n=l,2 COMPOUNDS 

Lets start by picturing the interacting electronic sys- 
tems within Sr2Ru04 and Sr3Ru207 at equilibrium for 
H=0. Besides the different number of Ru ions in the 
primitive unit cell, there are important differences in the 
crystal symmetry (see Fig. [T|). The n=l compound has 
ideal tetragonal symmetry (space group lA/mmm) with 
the fourfold rotation C^ around the c-axis. On the con- 
trary, the n=2 ruthenate shows orthorhombic symmetry 
(space group Bbcb), whereby the RuOg octahedra dis- 
play a small rotation of 6.8° in the a&-plane^, resulting 
only in the twofold symmetry element C2- Basic struc- 
tural and electronic differences may also be understood 
from an V2x V2 reconstruction of the Sr2Ru04 unit cell 
within the a5-plane. 

Figure [5] depicts the band structure and the local 
Ru(4(i) density of states (DOS) of the two layered com- 
pounds as retrieved from conventional LDA calculations 
without SOC. The low-energy regime is dominated by the 
t2g manifold of the Ru(4(i) shell and can be downfolded 
to associated Wannier-like states. In both systems the 




FIG. 3. (Color online) Sr2Ru04 quasiparticle band structure 
within standard LDA as well as with including SOC and cor- 
relations. 

Ir 





"•r M X 


r 


Z 


N M r 


1.5 


— "i,, 


i 


i 


- 


1.0 


— - d_. 


fi 
' 1 i 


li 


- 


0.5 




m 


|\ 


0.0_ 


l^w,J^V^ 


> ...iv 


i<^ 


4 -3 -2 


-1 n 

E-Ep (eV 


1 


2 1 




FIG. 2. (Color online) Comparison of band structure (top) 
and local Ru(4d) DOS (bottom) on the LDA level between 
Sr2Ru04 (left) and Sr3Ru207 (right) when neglecting SOC. 
The red curves in the top pictures display the Wannier-like 
dispersion as obtained from the MLWF scheme. The angular- 
momentum resolution is obtained from projecting the Bloch 
states onto cubic harmonics with a radial extension of 2.0 a.u.. 



FIG. 4. (Color online) Same as Fig. [S] here for Sr3Ru207. 
The lower panel provides a blow up of the dispersions close 
to the Fermi level. 





Jh=0.20 


Jh=0.35 


Jh=0.50 


U=1.5 
U=2.0 
[7=2.5 


-13.9 
-13.6 
-11.1 


-7.8 
-6.5 
-3.8 


-4.4 
-2.8 
-0.6 



TABLE I. Energy (in meV) of the highest occupied band at 
the X point in Sr3Ru207 for different U and Jh combinations 
(both in eV). In LDA without spin-orbit coupling that energy 
amounts to -22.8 meV. 



latter leak into the oxygen-dominated block of bands. 
The overall t2g bandwidth of the bilayer compound is 
slightly smaller (^3.1 eV) than for the single-layer sys- 
tem {^ 3.4 eV). The eg contribution close to the Fermi 
level is minor for Sr2Ru04, however the d^2_y2 charac- 
ter is non-negligible close to Ef for Sr3Ru207fi^ as may 
be seen from the local DOS. Note that this contribution 
is included in the effective t2g-MIJWF construction J^ In 
this respect it is important to record that we use in the 
following the (xy, xz, yz) terminology, though the cor- 




FIG. 5. (Color online) Interacting Fermi surface for Sr3Ru207 
in the basal ^^=0 plane from extended-LDA including SOC 
and electronic correlations. The black square marks the BZ 
cut and the labeling of the different sheets is according to 
Ref . fii. 



responding orbitals are only f2g-like in the sense of the 
present minimal MLWF construction. Note that these 
orbital functions (as true low-energy states) are point- 
ing inbetween the Ru-O-Ru bonds of the in-plane square 
lattice (compare also Fig.IHl). 

In the single-layer compound the energetics of the 
t2g manifold is split into {dxz, dyz) and dxy The 
dxz, dyz orbitals are truly degenerate, with crystal-field 



splitting Axy=ex 



= 113 meV to the d^y Wan- 



nier level. With very small quantitative differences, 
Sr3Ru207 has quasi-degenerate dxz, dyz levels (A=0.3 
meV) and A2;j,=115 meV. Hence it has to be empha- 
sized that there is already a small but nonzero splitting 
between dxz and dyz, giving rise to nominally slightly 
different local occupations. Note also that the strong 
local-DOS differentiation between [xz, yz) and xy does 
no longer hold in the bilayer compound. Albeit the dxy 
Wannier level is always lower in energy, the local orbital 
electron fillings in LDA for {dxy, dxz, dyz) are (1.24, 1.38, 
1.38) in the case of n=l and (1.40, 1.30, 1.30) for n=2. 
Thus there is a change in the t2g occupation hierachy 
between both ruthenates due to band-dispersion effects. 
Right at the Fermi level, the DOS of Sr2Ru04 is close 
to a van-Hove singularity slightly above ep, while in the 
case of the bilayer the Fermi level is located in a large dip 
of a complicated multi- valley-peak DOS at low-energy. 

Solving the problem posed by the minimal Hamil- 
tonian ([1]) for each system results in modifications in 
the low-energy dispersions. Figure [3] shows the t2g- 
like QP band structure of Sr2Ru04 from the extended- 
LDA treatment. One observes the expected combined 
main features already known from the existing separate 



SOG2^,30^ and correlated^SiSai^dl studies. Namely the 
lifting of degeneracies, e.g. close to the F point for the 
bands with dominant xz,yz character, also resulting in 
now avoided band crossings, e.g. close to the X point, 
when including SOC. Electronic correlations lead in the 
present approximation to band-narrowing and -shifting. 
Lifetime effects as well as incoherent spectral weight can 
not be retrieved within RISE at saddle-point. But the 
method captures very well the slightly modified Fermi- 
level band crossings in e.g. the TM direction as well 
as the shift of the van-Hove singularity towards the M 
point J^k 



lisa The band renormalizations are substantial, 
however not quite as strong as obtained within DMFT 
calculations with more elaborate frequency dependence 
of the self-energy2Si22ii^ to match the ARPES measure- 
ments. In the present modeling quantum fluctuations are 
missing and to obtain a ratio T?T,*/?TiLDA of the order of 3- 
4 for Sr2Ru04, in good comparison with photoemissioiii^ 
and dHvAi^ data, a value C/>3 eV would be needed. 

Accordingly, Fig. H] depicts the comparison between 
the quasiparticle bands (including SOC and correlations) 
with the conventional LDA dispersions for the bilayer 
Sr3Ru207. Due to the reduced symmetry the level of 
complexity is now surely raised. The complicated low- 
energy manifold with its rather flat bands already on the 
LDA level now shows significant effective-mass renormal- 
ization and additional splittings within the extended elec- 
tronic structure examination. From experiment , ^^i^-^ the 
ratio m*/TOLDA is effectively on the order of 6, thus even 
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FIG. 6. (Color online) Individual weights of the t2g Wannier 
orbitals on the low-energy QP bands from the extended-LDA 
calculation. 
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FIG. 7. (Color online) Interacting quasiparticle DOS for 
Sr3Ru207 compared to the LDA result. The QP content in- 
cluding SOC and correlations is scaled with Z=Q.7 to account 
for the modified integrated spectral weight resulting from the 
slave-boson framework. The bottom panel shows a blow up 
around the Fermi energy. 



larger than in the single-layer compound. Most interest- 
ingly, the blow up of the computed Fermi-level neigbor- 
hood renders it obvious that especially close to the X 
point severe changes take place. For instance a two-fold 
band approaching from M and splitting when passing 
X is strongly shifted to a low-energy of about 7 meV in 
the occupied part. Importantly, the latter small scale 
may only be reached by the combination of SOC and 
explicit Coulomb interactions. While the spin-orbit cou- 
pling is responsible for the appearance of split-up bands, 
the electronic correlations provide the narrowing and fur- 
ther shifting towards ep ■ Increasing U and Jh leads to an 
even closer placement nearby the Fermi level (see Tab.|I|. 
Hence the present approach is truly capable of describ- 
ing the very low-energy scale the bilayer compound is 
famously known for. Furthermore the key features of the 
named dispersive structure around X with its local min- 
ima and maxima are in close agreement with results from 
ARPES studies by Tamai et al.J^ The Fermi- level cross- 
ings in this region of the Brillouin zone (BZ) give rise to 
the so-called 72 pocket, which is under strong suspicion 
to play a vital role in the peculiar metamagnetic behav- 
ior of this compound4^ The complete Fermi surface (FS) 
obtained from our extended-LDA calculations exhibited 
in Fig. [S] is in good accordance with the one determined 
from photoemission^ Although however from our study 
the C|-symmetry character seems vital in the fermiol- 
ogy, but note that the experimental data in Ref. [l6| is 
symmetrized along the TX direction. 

An important question concerns the band characters at 
low energy in order to touch base with the local-orbital 
viewpoint. Figure [S] therefore shows the so-called fat- 
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FIG. 8. (Color online) Histogram showing the occupation 
probabilities of calculated multiplets in Sr3Ru207 ordered by 
the total angular momentum J. 



bands (orbital weight proportional to an artifical band 
broadening) for the t2g manifold. Its easily seen that 
the bilayer system is far from being a textbook exam- 
ple when it comes to attributing bands to a certain az- 
imuthal quantum numbers, since e.g. the supposingly 
relevant bands close to X are of strong mixed t2g char- 
acter. Nonetheless small asymmetries may be identified. 
The topmost occupied band at X has somewhat more 
dxz than dyz weight, true also for the lowest one in the 
given energy window. In between there is only one with 
more dyz character. The d^y orbital has notably overall 
the same order of weight in this region of the BZ as its 
out-of-plane companions. Interestingly, from the strong 
weight along TA the dxy character seems to dominate the 
propagation perpendicular to the RuOg planes. 

The low-energy scales may also be confirmed from the 
quasiparticle DOS plotted in Fig. [71 Therefrom its again 
obvious that the states close to ep are strongly pro- 
nounced in extended-LDA, shifting prominently to the 
Fermi energy. The LDA-DOS exhibits a smaller peak at 
the Fermi level within a valley of ~60 meV width. With 
the additional interactions that feature is strengthened 
and importantly a peak in the low-energy occupied region 
is sharpened and shifted towards ep, being located at ~10 
meV. These findings of increased spectral weight below ep 
within a meV range is in accordance with photoemission 
studiesi^ and also supported from specific-heat dataj^ 

Finally, Fig. [5] depicts the occupation probabilities 
of the local t2g-based multiplets according to the con- 
verged slave-boson amplitudes of the lattice calculation 
in the metallic state. There are various sizeable mul- 
tiplet weights, non-surprisingly with an overall domi- 
nation of the ones from the four-particle sector. The 
atomic ground-state multiplet with L=S=J=l, an or- 
bital and spin triplet, has also the largest weight in the 
itinerant regime. Note that the deviations of J from 
the ideal values results from the small Cg weight inter- 
mixing within the LDA-derived Kohn-Sham Hamiltonian 



based on the low-energy t2g-\ike orbitals. For Sr2Ru04 
the corresponding picture looks very similar, also there 
the L=S=J=1 multiplet has the maximum weight. 



IV. BILAYER RUTHENATE IN APPLIED 
MAGNETIC FIELD 



Since our approach is in the position to correctly ad- 
dress the equilibrium low-energy correlated electronic 
structure, we expect a qualitatively meaningful descrip- 
tion of Sr3Ru2 07 with applied magnetic field H. The 
Zeeman-type local interaction is now included together 
with the spin-orbit and Coloumb interactions. All are 
adequately treated due to the generality of our formal- 
ism, also allowing for arbitrary field directions. Note that 
the crystal structure is constructed such that the in-plane 
square lattice evolves along the Ru-O-Ru bonds, while 
the original x-, y-axes point inbetween these bonding di- 
rections (see Fig. [9]). In the present investigation the 
field direction is modified within the xz plane, i.e., the 
tilting of H takes place along the diagonal of the square 
lattice which also agrees with the orthorhombic a axis. 
Notably the x direction in real space corresponds to the 
FM' (M rotated by 90° degrees) direction in reciprocal 
space and therewith the FX(X') direction amounts to 
QP propagation along the Ru-O-Ru bond on the rods of 
the square lattice. The resulting net magnetic moment 
M per Ru ion is composed of spin and orbital- momentum 
parts, i.e. Af=/iB(g J-H) and is here computed on the lo- 
cal level from the self-consistent slave-boson amplitudes 
(see eq. (HU) and not from fc- integrating the associated 
QP contributions. While in effective single-particle cal- 
culations both approaches would yield identical results, 
in slave-boson theory those numbers may in principle dif- 
fer since the QP part only carries the contribution to the 
itinerant character of the electron. Via eq. (jlip it is en- 
sured that the occupation number of the physical electron 
is retrieved. 





FIG. 9. (Color online) Left: view along the c axis of sketched 
Sr3Ru207 in order to clarify the invoked cartesian coordi- 
nate system with Ru ions in (blue/dark) and O ions in 
(red/grey). Right: applied magnetic- field direction where the 
z axis equals the crystallographic c axis. 



Because of the low-energy scales involved in this prob- 
lem, especially for the close-to-realistic investigation of 
the magnetic behavior the numerics is however in any 
case far from simple. The complicated local interacting 
Hamiltonian asks for about 2000 slave-boson amplitudes 
to solve for, interlinked with a 24x24 QP- Hamiltonian 
problem (four Ru ions with three orbitals per ion, allow- 
ing for spin degrees of freedom within each orbital) on a 
properly dense fc-point mesh. 



A. The case H | c 

We first choose the crystallographic c-axis (i.e., the 
z direction) perpendicular to the Ru02 planes for the 
magnetic-field direction, i.e. H=_ff c. Figure [10] shows 
the evolution of M with increasing field strength together 
with the obtained total free energy for the bilayer system. 
Both curves display a rather non-trivial behavior. In the 
following the values of H are given in meV. Up to H=5 
the value of M rises linearly (region I), with however neg- 
ative free-energy curvature, hinting towards unfavorable 
field penetration. In the range 5<H<23 a stronger rise of 
M occurs, followed by a non-monotonic behavior (region 
II). In that second region the free energy is first deplet- 
ing and then again rising along with the non-monotonic 
part, overall evolving with the positive curvature of a sta- 
ble phase. Lets also note that it appears as if in the first 
part 5<H<10 the M evolution as well as the free-energy 
curve display some non-trivial modulation. Finally for 
even larger H (region III) the magnetic moment enhances 
further with close-to-linear development. In III the free 
energy follows a novel parabolic shape at higher values 
than in II. The described evolution reveals the metam- 
agnetic behavior of M, with strong resemblance to the 
experimental data^Sii^. Notably, the observed metamag- 
netism is again only obtained in the present computations 
if both, SOC and electronic correlations, are included. 
Due to the strong signatures in the free energy, the tran- 
sition between regions I/II as well as regions II/III are 
clearly of first order. We however did not investigate the 
transition orders in a more elaborate fashion (e.g. via 
computing the Hesse matrix). 

To connect these global results to the orbital degrees of 
freedom. Fig. [TT] depicts the orbital-resolved local occu- 
pations and contributions to the magnetic moment. With 
increasing field the dxy filling shrinks, while the one for 
dxz, dyz grows. Thus notably there is an inter-orbital 
charge transfer from dxy to dxz, dyz with magnetic field. 
Whereas for H=Q a marginal filling difference between 
the latter orbitals is observed, with growing magnetic 
field the occupations of the quasi-degenerate levels more 
or less align. The respective orbital contributions to the 
magnetic moment within the t2g manifold are strongly 
varying. While the dominant dxy part shows substan- 
tial paramagnetic response with field, the generally much 
smaller dxz, dyz terms exhibit intricate behavior. They 
start with flat, nearly constant minor diamagnetic re- 
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FIG. 10. (Color online) Sr3Ru207 in applied magnetic field 
(measured in meV) along z. Left: free energy vs. net magnetic 
moment, right: net magnetic moment per Ru ion. 
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FIG. 11. (Color online) Orbital-resolved local occupation 
(top) and contributions to the Ru net magnetic moment (bot- 
tom) with magnetic field. The inset shows a blow-up of the 
dxz, dyz curves. 



sponse for small magnetic field and both only turn into 
weak paramagnetic characteristics at H^li. Interest- 
ingly, this PM behavior shows small differences in the 
amplitude for dxz and dyz (with Mxz>Myz), with even 
an observable sudden increase in that difference when the 
II/III transition occurs. As the sole d^y response shows 
no definite MM signals, the para/dia discrimination in 



the orbital response of dxz ■, dyz appears as a key micro- 
scopic building block for the MM behavior. 

The orbital discrimination in the magnetic behavior al- 
ready provides a first clue to the MM puzzle of the bilayer 
system. A second important insight originates from fol- 
lowing the development of the QP bands with magnetic 
field, presented in Fig. [T^l It is seen that increasing H 
amounts to intriguing changes in the low-energy mani- 
fold, most notably to relevant shifts in the peculiar QP 
structure close the X point. For H=4, non-surprisingly 
there are further band splittings along X compared to 
the case of zero magnetic field and the lowest-energy p- 
shaped band at this fc-point is now placed just below 
s-p. A discrimination between the 72 pockets around X 
and X' is furthermore clearly visible. Going to H=6 the 
named band is locked to the Fermi level and 72 at X 
has opened towards the BZ boundary where as 72 at X' 
has shrunk. For H=10 the former band is above ep and 
both 72 pockets are opened. At H=18 the second- lowest 
band at AT, X' crosses the Fermi level and a minor pocket 
structure reappears at X' . Morover the a2 sheet starts 
to become increasingly distorted along FA (A') and also 
begins effectively shrinking with growing field. The latter 
signatures are strengthened for H=26 with the additional 
onset of hybridization between ai and 6. The pockets 
close to A, A' remain both opened in that large-field re- 
gion III, with now three low-energy bands having crossed 
ep at A, A'. Along with these changes, the QP-DOS 
of course runs through several peaks at the Fermi en- 
ergy, but evidently exhibiting an evolution different from 
a pure shifting of the H=0 structure. Thus from the 
Bloch perspective the picture of Lifshitz transitions un- 
derlying the magnetic response emerges. Various authors 
have already pointed out the importance of van- Hove sin- 
gularities crossing the Fermi level and we here can verify 
this mechansim based on the complete realistic starting 
point. The fact that not only the 72 sheet but also the 
inner sheets, most notably a2, may play a vital role in 
the MM response was also retrieved in a recent experi- 
mental study employing spectroscopic imaging scanning 
tunnelling microscopy<^ 

The applied magnetic field leaves also some signa- 
tures in the orbital-dependent electronic self-energy (see 
Fig. [Tn)l . All orbital sectors display the expected splitting 
in the QP weight Z and static self-energy S^*^' due to the 
spin-filling inbalance with larger magnetic polarization. 
The splitting for dxz, dyz is weaker and especially minor 
at small field where the nearly constant diamagnetic re- 
sponse occurs. Overall there is no strong modification 
of the correlation strength with H, the MM signatures 
show up somewhat stronger in the QP weight that is 
associated with band renormalizations. From the calcu- 
lation, the value for Zxy is slightly larger than for the 
remaining two <2g orbitals. But because of the intriguing 
hybridizations (compare Fig. [6]) no trivial relation may 
be drawn therefrom in view of the respective renormal- 
ized effective inasses on the various Fermi sheets. Our 
RISB formalism is furthermore able to reveal the local- 
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FIG. 12. (Color online) Development of the QP band structure, DOS and Fermi surface with magnetic field along the z 
direction. The (blue/grey) circle on the right panel always marks the respective H value. 



multiplet behavior of the correlated system. Via the 
slave-boson amplitudes the method allows to evaluate the 
occupation probability of a given eigenstate of the local 
Hamiltonian ([2]) within the complete itinerant solution. 
For instance, Fig. [T3] depicts the splitting characteristics 
of the *S'3/2 multiplet from the three-particle sector as- 
soicated with orbital momentum I/=0, i.e. having J as 
purely spin-defined. As exptected, the J^ degeneracy is 
lifted for H^Q, showing non-trivial signature close to the 
phase transitions between the different region I-III. Sure 
enough, the four-particle sector is most dominantly oc- 
cupied for the Ru(4(i) t2g shell, but the multiplets there 
do not exhibit a very conclusive behavior with applied 
H. 



B. The case H |f c 

In principle our approach works for arbitrary polar an- 
gles 9 between H and the c-axis of the system. How- 
ever since the computations are rather expensive we have 
chosen here only two additional specific values (besides 
6'=0), namely 6'=20° and 6'=90°. The azimuthal angle is 
put to zero, i.e. the latter in-plane magnetic field points 
along the x direction inbetween the Ru-O-Ru bond of 



the square lattice (see Fig. [S]) . The evolution of the net 
magnetic moment per Ru ion for the two new field angles 
together with the respective free-energy plot is displayed 
in Fig. 1151 When directly comparing the magnetizations 
for 9=2Q° with the former 6'=0, one first realizes that for 
a given H the value of M is increased. This seems to be 
in line with the experimental data from the work of Grig- 
era efc aii^ showing an enhancement of the real part of 
the differential susceptibility at the MM transition with 
9 (note that is defined as the angle between field and 
the afe-plane in Ref. [221 ) . The overall phenomenology of 
M{H) for 0=20° is still rather similar to the case H||c. 
Note that although the upper first-order transition hap- 
pens at larger magnetic moment for 0=20°, the value 
corresponds to nearly the same magnetic-field strength 
H . On the contrary, for H_Lc along x the overall char- 
acteristic is qualitatively different. After a near linear 
rise for H>5, close to H=16 a sudden change of slope 
for M{H) takes place (with a possible signature in the 
free energy) and the magnetic moment continues again 
nearly linearly. The upper first-order transition is not vis- 
ible anymore in the free-energy curve. This observation 
shows clear resemblance to the experimental findings of 
singular behavior for in-plane magnetic field compared to 
strong out-of-plane Hi^ However from the present com- 
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FIG. 13. (Color online) Orbital-dependent QP weight (right) 
and static self-energy (left) with H . 
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FIG. 14. (Color online) Splitting of the ''53/2 inultiplet in the 
applied field. 



putations we can not draw a definite conclusion concern- 
ing the shift of the phase boundaries with respect to H 
and Q. We are nevertheless in the position to shed more 
light on the angular-dependent differences by showing in 
Fig. [16] the orbital contributions to the occupation and 
the magnetic moment for 0=20°, 90°. For the smaller 
Q the orbital-resolved differences compared to H along c 
are marginal, with mainly an obvious occupation differ- 
ence between dxz-, dyz (with nxz>nyz) for H>22 which 
was absent before. Comparing the results for H along 
X to the case along z shows again clear signature. The 
stronger occupation of d^z occurs at already small mag- 
netic field and becomes substantial at large H. Further- 
more the former area of diamagnetic response from d^z , 



dyz has nearly vanished, both orbital responses are nearly 
indistinguishable from zero within the accuracy for small 
H. Yet at least for H>15 the PM response with differ- 
ent amplitude is clearly observable. Thus the dia/para 
competition in dxz , dyz that seemed to be crucial for the 
MM behavior for H along z is nearly absent, explaining 
the qualitative difference between in- and out-of-plane 
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FIG. 15. (Color online) Free energies (left) and net mag- 
netic moment (right) for ^=20° (top) and 6=90° (bottom), 
compared to the moment for H along the z direction (6=0°), 
respectively 




FIG. 16. (Color online) Orbital-resolved occupation and con- 
tributions to the magnetic moment as in Fig. Illl here for 
61=20° (left) and 6=90° (right). 
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FIG. 17. (Color online) As Fig. [H here for 61=20° (left) and 61=90° (right). 



field. Coming back to the phase-region shift with field 
and angle, and fixing such a shift to the dia/para crossing 
for dxz , dyz ; one might observe that this crossing indeed 
shifts to the left, i.e. towards smaller H . However again 
the resolution is not accurate enough to render a unique 
statement concerning that question. The impact of the 
finite angle 9 shows also up in the changes of the QP 
states with iJ, as documented in Fig. [TT] While again 
for 0=20° there are no major qualitative differences com- 
pared to the case of H along c (e.g. also here the three 
low- lying bands at X cross the Fermi level with field), the 
pure in-plane field leads to clear modifications. Namely, 
the two originally nearly degenerate lowest-energy bands 
at X do not split with H, but cross e-p together. In line 
with this, the four 72 pockets in the BZ behave coherently 
for all field strengths and also the 5-ai hybridization at 
large H occurs now in a fourfold manner. Importantly 



the third low-energy band (now also still degenerate with 
the fourth one) remains below the Fermi level within the 
range of the studied magnetic field. Thus H along the x- 
axis leads to an avoided lifting of degeneracies, resulting 
in qualitative different magnetic response. However note 
that this on the other hand does not imply that the sym- 
metry between dxz , dyz is enforced, since the local orbital 
discrepancy is strongly increased for the sole in-plane 
field (compare Fig. [T6| . Such selection-rule constraints 
depending on the magnetic-field direction show up also 
prominently when it comes to local-multiplet splittings, 
as shown in Fig. [TH] for our example of the L=0 ^83/2 
many-body state from the three-particle sector. While 
the Jz splitting for 0=20° remains vital, in the case of 
H along X the splitting is mostly absent. For further 
clarification. Figure [11] depicts the resulting angle be- 
tween applied field and net magnetization, rendering it 
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FIG. 19. (Color online) Angle between applied magnetic field 
H and resulting net magnetic Ru moment for 0=0, 20°, 90°. 



clear that only at large enough field strength the moment 
aligns along H . In this respect it appears as if for H_Lc 
the moment is somewhat more easily forced into the field 
direction. 



V. SUMMARY AND DISCUSSION 

The physical content of this work is twofold. First 
the interplay of spin-orbit coupling and local Coulomb 
interactions was studied at equilibrium for the n=l,2 
layered strontium ruthenium oxides belonging to the 
Ruddlesden-Popper family. Thereby we started from the 
realistic low-energy Kohn-Sham dispersion as obtained 
from Wannier-downfolding the bands from state-of-the- 
art LDA calculations. Then notably the multi-orbital 
many-body effects (and its interlinking with spin-orbit 
effects) were treated beyond simple Hartree-Fock mean- 
field by utilizing proper self-consistent renormalizations 
due to strong correlations within rotationally invariant 
slave-boson theory at saddle-point. 

For both compounds it became evident that SOC and 
strong correlations together are important to account for 
the detailed low-energy electronic structure at small tem- 
peratures. Already in Sr2Ru04 the local Coulomb cor- 



relations effectively renormalized the spin-orbit interac- 
tion, leading to enhanced band-splittings for the t2g man- 
ifold close to the Fermi level. Its is therefore surely ex- 
pected that the unconventional superconductivity at low 
T has to be addressed by treating these both interac- 
tion types on equal footing.— For the bilayer compound 
Sr3Ru207 such an approach was shown to be essential in 
order to describe the intriguing low-energy quasiparticle 
band structure and density of states in close resemblance 
to existing ARPES and specific-heat measurements.— 1^ 
The renormalized spin-orbit split bands give rise to ex- 
tremely small energy scales, whereby close to the X point 
in the BZ an especially rich structure appears. However a 
straightforward decomposition of the complicated renor- 
malized band structure into distinct dxy-, dxz-, or d^^-like 
bands seems difficult, the system looks like an intricate 
multi-orbital system where subtle differences in the or- 
bital contributions eventually play a crucial role. 

Beyond the equilibrium study, a straightforward exam- 
ination of the bilayer ruthenate in applied magnetic field, 
based on the complete multi-orbital t2g Hamiltonian in- 
cluding the Zeeman term in the presence of SOC, was 
presented. Depending on the magnetic-field direction, 
metamagnetic transitions in line with first-order phase 
transitions were verified. Taking a local viewpoint in that 
metallic system, the competition between paramagnetic 
(dxy) and diamagnetic {dxz, dyz) contributions appears 
to play a significant role for the MM phenomena. More- 
over an orbital charge-transfer from dxy to dxz, dyz with 
increasing H was observed, though the orbital fillings 
between the latter two t2g orbitals only seem to deviate 
with larger field angle 9. Concerning the itinerant QP 
states it became evident that Lifshitz transitions close to 
the X point (i.e. around the 72 pocket) may partly be 
blamed for changes in the free energy and accompanying 
Fermi-surface reconstructions across the MM transitions. 
But moreover a substantial change of the a2 sheet and 
an additional ai-S hybridization is seen with larger mag- 
netic field. While the region around X corresponds to 
propagation along the Ru-O-Ru bond with strong short- 
range variation, the latter inner sheets account for rather 
isotropic in-plane transport with only long-range varia- 
tion in the respective QP wave function. Thus the physics 
of the MM transitions involves directional, short-range 
processes at smaller H and incorporates long-range mech- 
anisms at larger H . In addition, symmetry changes seem 
to take place between certain regions of the BZ. Albeit 
the overall C2 symmetry of the total FS remains sta- 
ble with H, the original C|-like symmetry between the 
72 pockets is disturbed towards C2 in the central MM 
phase region II, if H||c holds. However that symmetry- 
change is absent for H_Lc and a well-defined bounded 
MM region can not be identified. Thus the present cal- 
culations reveal the qualitative differences between out- 
of-plane and in-plane magnetic-field in accordance with 
experiment (see Ref. |Tl| for a recent review). 

There remain however open questions. For instance, 
we may not draw definite conclusions on the shifting of 
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the MM transitions with field angle, whereas in some 
model studies^i^ the shift of the MM phase region to 
lower magnetic fields with 9 in line with experimental 
work was verified. Besides several other possible reasons, 
this angle-dependent behavior might also be sensitive to 
the specific choice for the magnitude of the interaction 
parameters U , Jh and A within the local Hamiltonian, 
also in conjunction with the nesting properties between 
the various spin-polarized Fermi sheets 4^ 

Another important point concerns the appearance 
of nematic order, revealed in transport studies for 
Sr3Ru207 to escort the MM region)^ We evidently see 
symmetry changes in the FS geometry and standard 
representations of nematic order parameters for a cer- 
tain angular-momentum ordering channel I of the form 
Mi=22-^^n{k) cyi-p[ilip{]^)] (see e.g. Ref. 1671 for a review) 
display non-trivial behavior with H depending on I and 
on the number of included bands. However the com- 
puted data does not exhibit convincing evidence for a 
well-defined quantification of nematicity along that def- 
inition. Note that we also did not incorporate such an 
symmetry-breaking (i.e. forward-scattering) term explic- 
itly in the Hamiltonian, as done in some model stud- 
jgg^^j4o^j44 ^^^ j^ might therefore be possible that 
we miss additional (energetically favorable) symmetry- 
breakings in our mean-field approach. For the present 
problem, the nematic order parameter has frequently also 
been defined via the filling difference between the d^z , dyz 
orbitals42iiL Yet this definiton seems dangerous, since 
the C2 symmetry of Sr3Ru207 is broken already by the 
equilibrium crystal structure and LDA calculations re- 
veal such a nominal filling difference for H=0. On the 
other hand our extended-LDA calculations point towards 
an initial alignment of these sub-orbital fillings with mag- 
netic field. Only for larger {H, 9) a true d^z, dyz filling 
differentiation occurs. In this respect its also noteworthy 
that here the magnetic-field tilting towards the a6-plane 
takes place inbetween the Ru-O-Ru bond, i.e. along the 
diagonal a-axis of the in-plane square lattice. It would 
therefore be interesting to check further additional in- 
plane directions, especially along Ru-O-Ru. 

This brings us to possible extensions of the current 
work on Sr3Ru207. For present numerical reasons, the 



different Ru ions in the bilayer unit cell were assumed 
equivalent by symmetry (as true for the equilibrium crys- 
tal structure). But we easily expect that the MM phase 
regions would generally benefit from such a symmetry 
breaking within the unit cell. If orbital-liquid scenarios 
or the physics of domain structuresii*^ in the compound 
are vital, one has to come up with even more sophisti- 
cated real-space picturings. Inter-site Coulomb interac- 
tions are here neglected, but may also be a source for the 
observed symmetry breakings 4^ Furthermore since this 
work was performed by using a postprocessing scheme 
to existing LDA calculations, elaborating on a complete 
charge-selfconsistent approach with the proper feedback 
of the electronic self-energy onto the Kohn-Sham charge 
density could surely enhance the MM response. Ac- 
counting for finite-temperature effects is an additional 
further important aspect in order to reveal the intri- 
cate thermodynamics of the MM region<^ Last but not 
least, at the moment the method relies on mean-field 
theory. Since in other strongly correlated materials it is 
already found that intricate self-energy effects especially 
take place close to van-Hove singularities^^ local quan- 
tum fluctuations should be included in future Sr3Ru207 
studies. Non-local quantum spin fluctuations may have 
relevant impact on the low-energy physics, since the ma- 
terial is prone to magnetic order. Nonetheless the present 
realistic formalism yields promising results and shows 
that extended-LDA calculations are in principle capable 
of addressing the challenging low-energy physics of the 
layered ruthenates. 
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